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ABSTRACT 

As we know, let alone to find the determinant of infinite matrix, it is difficult to find the determinant of some n x 
n matrixes by the usual methods like, the cofactor method and Crammer’s rule. But now we will show how to find the 
determinant of some n x n matrices and how to find the determinant of some infinite matrix by using Exponential 
Generating Functions. In this paper we will consider matrices having 1,2, 3, 4, 5... on the supper diagonal, 0's on the upper 
and identical entries on each diagonal below the supper diagonal. Here we will try how to obtain the determinant of n x n 
upper left corner sub matrix of a given infinite matrix by introducing Exponential Generating functions of some sequences 
and how to get a sequence by calculating the determinant of n x n upper left corner sub matrix of infinite matrix, we will 
also check the correctness of the determinant by using Numerical method. 

KEYWORDS: Infinite Matrix; Determinant of Matrices; Exponential Generating Functions; Sequences; Sub Matrix 

1. INTRODUCTION 

To understand the whole work, it is better to know about a matrix, determinants. Generating functions and some 
sequences. So we will discuss these terms before the actual work. 

What are Generating Functions? 

One of the main tasks in combinatory is to develop tools for counting. Perhaps, one of the most powerful tools 
frequently used in counting is the notion of Generating functions. Generating functions are used to represent sequences 
efficiently by coding the terms of a sequence as coefficients of powers of a variable x in a formal power series. Generating 
functions can be used to solve many types of counting problems, such as the number of ways to select or distribute objects 
of different kinds, subject to a variety of constraints, and the number of ways to make change for a dollar using coins of 
different denominations( Discrete Mathematics and Its Applications, Seventh Edition, Kenneth H. Rosen Monmouth 
University (and formerly AT&T Laboratories page 537). 

A generating function is a continuous function associated with a given sequence. For this reason, generating 
functions are very useful in analyzing discrete problems involving sequences of numbers or sequences of functions. 
(Generating Functions and Their Applications Agustinus Peter Sahanggamu MIT Mathematics Department Class of 2007 
18.104 Term Paper Fall 2006 page 1) 

In mathematics a Generating function is a formal power series whose coefficients encode information about a 
sequence { a n } that is indexed by the natural number n. Generating functions can be used to solve determinants of some nxn 
and then an infinite matrix by relating the terms of the sequence for which we get a generating function to the determinant 


Impact Factor(JCC): 2.7341 - This article can be downloaded from www.impactjournaIs.us 






44 


Masreshaw Walle Abate 


of an upper left comer nxn matrix of an infinite matrix.. Even though there are various types of Generating functions, 
in this paper, we introduce the idea of Exponential generating Functions. 

Definition 

The exponential generating function A(x) of a sequence { a n } 00 n is defined as 


2 3 

XXX 


Z -/v -/v -/v -A- 

a r = a Q + a { — I - a 2 1- a 3 1- .... a r — .... the exponential generating functions have 

,. =0 r! 1! 2! 3! r! 


the following property. 

Lemma 


If A(x) is the exponential generating function of a sequence \ci n } j|=Q then 


2 n= 0 a„ + i— = A ( x ) 


Proof: Differentiate both sides of the definition of A(x) in 


X 

r= 0 n\ 


2 3 

XXX 


Z A -A. -A- -A A 

a„ — = a n +a,— + a 0 b a, (-.... a ,. — .... 

" u 1 1! 2 2! 3! r! 


A'(x) = mflm — =Sm=! “m =S=0 “n+1 ^T 

Some Examples of Generating Functions of Some Sequences 

i) (1,1,1,...) <— > 1 + x + X 2 + X 3 + ... = x" = is the ordinary Generating function 

n-0 1 “ X 

.111 . X x XX 

ii) (1,1, — , — , — > — — e =l + X-\ 1 1- .... is its exponential Generating function. 

2! 3! 4! ^ r! 2! 3! 


iii) ( 1, 2, 3, 4...) <H> l+2x+3x 2 +4x 3 +... : 


(1-x) 2 


is the generating function for counting numbers. 


2 3 2 2 3 3 

iv) The generating function for the sequence (1, k, k~, k ...), where k is an ordinary constant is 1+kx +k~x + k x +... = 

1 


1 - kx 


v) ( 1,-1, 1,-1...) ol-x+ x 2 — x 3 + ... = 

vi) ( 1, a, a 2 , a 3 ... ) 1+ ax + a 2 x 2 +...= 


1 


1 


1 - (-x) 1 + X 

1 


1 - ax 
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vii) ( 1, 0, 1, 0, 1... } <-> 1+ x 2 + x 4 + x 6 +... = 1 

1-x 2 

1.2 OPERATIONS ON GENERATING FUNCTIONS 

Let A(x) = a 0 +ai x+a 2 x 2 + a 3 x 3 +... and B(x) = b 0 +b j x+b 2 x 2 + b 3 x 3 +... be the generating functions for the 
sequences (a r ) and (b r ) respectively, then 

1.2.1 Product Rule (Convolution) 

Although multiplying two functions is as natural as adding them, as it is seen below, general term of the resulting 
sequence of multiplication is not the product of i th terms of the sequences. In fact: - 

A(x)xB(x) = a 0 b 0 + (a 0 bi+ a!b 0 )x +( a 0 b 2 + a 3 bi +a 2 b 0 )x 2 ....is the generating function for the sequence (c r ), where 

C|-— a 0 b 0 + aib r _i+ — a r _i bj+ a r bo, r = 0,1, 2, 3,... 

Proof 


To evaluate the product a(x) x B(x) let us use the following table. 

Table 1: Product Table 



If we follow the arrow, we get the required product 

E.g. If A(x) <-> < 1, 2,2,2, } ^- l — ' - and B(x) <-> ( 1,1, 1,...) <-> — — then 

1-x 1-x 

A(x). B(x) = ( — )( — — ) = l + X ^ < 1, 3,5,7, 9.... ) 

1-x 1-x (l - x)“ 

Using the product rule we have the following: 

(1-x) A(x) is the generating function for the sequence (c r ) where 
c 0 = a 0 and c,= a r - a r _ for all r >1 and 

A 

— ^ is the generating function for the sequence (c r ) where 
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c r = a 0 +a,+ a 2 +... + a r for all r. 


Remark 


2 3 n 

since e' = 1 + X H b — — 1- and 

2! 3! n\ 


r 2 r 3 r" 

e~ x = l-x + - — + ... + (-1)" — + . 

2! 3! n! 


x , — x 2 4 

C \ C X X 

Then we have = lH 1 b... = COS* and 

2 2! 4! 

e x ~e~ x A 3 A 5 

= x-\ 1 b... = Sin A 

2 3! 5! 

1.2.2 Addition Rule 

A(x)+B(x) is the generating function for the sequence (c r ) where c,= a,+ b r , r = 0,1, 2, 3... 

Proof 

A(x)+B(x) = (a 0 +a 1 x+a 2 x - + a 2 x +...) + (bo+bjX+b 2 x - + b 2 x +...) 

= ( a 0 + b 0 ) + (a; + bi)x+( a 2 + b 2 )x~ +(a 3 + b 3 .)x'+... 

= C()+ c^x + C 2 X - + C 3 X +... 

where c,= a r + b r , r = 0, 1,2,3... 

eg. if AO) = ^ = 1 + 2x + 4x 2 + 8x 3 + ••• ( 1,2, 4, 8, ... ) 

B(x ) = ^ = 3 = 3(1 + 3x + 9x 2 + 21 x 3 + ••• ) = 3 + 9x + 27x 2 + 81x 3 + ••• <-> ( 3,9,27, ... ) 

[ 3 4-9x 

then A(x) +B(x) = + = <-> ( 4, 11,31... ) 

1-2* 1-3* 1 - 5* + 6x 2 

1.2.3 Constant /Scaling/ Rule 

( ca 0 , ca b ca 2 , ca 3 ca n ... ) <->cA(x). 

Proof 

( cao, cai, ca 2 , ca 3 ... ) <— > ca 0 + caix+ca 2 x 2 +ca 3 x 3 +... 

=c(ao+a[X+a 2 x 2 +a 3 x 3 +... ) = cA(x) 


E.g. If < 1,2, 3,4,...) l+2x+3x 2 + 4x 3 +... : 


1 


( 1 -*)' 


and c=2, we have 
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— = 2 +4x + 6x 2 + 8x 3 +....^ (2,4,6,8,10,...) 

( l-x ) 2 v y 

1.2.4 THE DERIVATIVE RULE 

If ( a 0 a] a 2 ... ) <— > A(x), then ( a t 2a2, 3a 3 , ... ) <— >A'(x) 

Proof 

( ag ai, a 2 ) ^ ^ A(x) = ao + a^x + a 2 x~ + a 3 x +... 

=> A ( x ) = a[+2a 2 x + 3a 3 x 2 +... = A'(x) <-> ( ay2a 2 , 3a 3 2 ’... ) 

dx 

E.g. ( 1,1,1,.... ) <-> — — = 1 + x + x 2 +x 3 +...= A(x) 

l-x 


d ( l } 



dx 


U-^y 1 


= l+2x + 3x 2 +... VA ( 1, 2, 3, 4,... ) = A'(x) 

(i-0 


And x A'(x) VA ( 0, 1, 2, 3,... ) VA — 

(l-x ) 2 


Hence [xA'(x)]' = 


( V 

X 




+ x) 
(1+xY 


<— > ( 1,4,9,16,... ) (square number Sequence) 


Note:|Y"J Given a generating function A(x) we use|x n J A(x) to denote a n , the coefficient of x 
(270 Chapter 10 Ordinary Generating Functions) 

MATRIX AND DETERMINANTS 

Definition 


A matrix is a rectangular array of mn quanties ay 


i = l,2,3,....,»t & 


in m - rows and n- columns. It is called an m 


v j = l,2,3,...,/i 

Xn matrix or a matrix of order m X n and read as m by n matrix. The numbers ay are called the elements (constituents or 
coordinates or entries) of the matrix and we will denote the matrix by { ay } or A. The suffix ij of an element ay indicates 
that it occurs in the i th row and j th column, when n=m we call this a square matrix. For a square matrix nxn, if n-> co then 


the matrix is called an infinite matrix. 
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In Explicit form 



«n 

a 12 

a i3 

• • a u 


a 2\ 

a 22 

a 23 

■ ■ a 2n 

A = 

a 31 

a 32 

a 33 

■ ■ a 3n 



a n2 

«„3 

■ ■ a nn 


rows 


*- 




Columns 


NOTE 


1) If A — { ay } m x n and B— {by } n x p, and AB — C, then C — { Cy } mx p 


where Cy = ^ a tJ b 

k=l 


kj 


2) For any matrix 



MATRIX MUITIPLLCATION 

Two matrixes A and B are conformable for the product AB when the number of columns in A is equal to the 
number of rows in B. If A is an m x n matrix and B is an n x p matrix then their product AB is defined as m x p matrix 
whose (ij) th element is obtained by multiplying the element of the i th row of A in the corresponding elements of the j th 
column of B and summing the products so obtained. So the (ij) th element of the product AB is the inner product of the i th 
row of A and the j th column of B. 

DETERMINANTS 

Definition 

Determinant of a matrix A is a specific real number assigned to A It is denoted by det (A) or |A| Or for n > 1 the 

determinant of an n x n matrix A= (ay) along the first row is the sum of n-terms of the form + ay det Ay with plus and minus 
signs alternating where the entries a n ai 2 ai 3 ...a ln are from the first row of A. 
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In symbols, det A = an det An - ai 2 det A ]2 +... + (-1) n+1 a in det Ai„ = ^ (-1) l+ ’ a^ Ay and det Ay is the 

j = i 

determinant of the sub matrix which is obtained by removing the 1 st row and the j th column. 

Actually det Aij is called the minor of a^ and (-1) +J aij det Ajj is called the cofactor of ay. 

CRAMER'S RULE 

Let A be an invertible n x n matrix. For any b in R n the unique solution x of 

Ax = b has entries given by x , = w ^ ere ' = 1>2, 3,... n and A , (b) is the matrix obtained from A by 

replacing column i by the vector b. 

Proof 

Denote the column of A by a 1 ,a 2 a 3 ,... a n and the column of the n x n identity matrix I by e 1 ,e 2 ,e 3 ,... e n 
If Ax =b then the definition of matrix multiplication shows that 
A Ij(x) = A [ei,e 2 .e 3 ,... x... e n ] = [Aei,Ae 2 ... Ax... Ae n ] 

= [ai,a 2> ... b... a n ] = A ; (b) 

by the multiplicative property of determinants 

(det A) det 1 1 (x) = det A , (b) 

=> (detA) X j = det A , (b) 

_ x _ ctetAj(b) 

1 detA 

Determinant of a matrix can be obtained by the cofactor method or by using the Cramer's rule . But now we are 
Interested to show how to find the determinant of several matrices by using Generating functions. 

The matrices whose determinants we will be evaluating have all l's on the super diagonal . 0’s above the supper 
diagonal ; and identical entries on each diagonal below the supper diagonal, perhaps with the exception of the first column, 

2. DESCRIPTION OF THE METHOD 

In this topic we will see how to get a sequence from the given matrix by calculating the determinant of each 
upper left nxn square matrices of the given matrix. All matrices in this section will have 1, 2, 3, 4, 5, on the supper, 
0's above and identical entries on each diagonal below, perhaps with the exception of the first column. Hence we will 
equate the n th term of the sequence with determinant of each upper left nxn square matrices of the given matrix. We begin 
with a typical example as follows 
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Example: A 


10 

1 

0 

0 

0 

0 

0 . . . 3 

1 

0 

2 

0 

0 

0 

0 . . . 

1 

1 

0 

3 

0 

0 

0 . . . 

1 

1 

1 

0 

4 

0 

0 . . . 

1 

1 

1 

1 

0 

5 

0 . . . 

1 

1 

1 

1 

1 

0 

6 . . . 

1 

1 

1 

1 

1 

1 

0 . . . 

1 

1 

1 

1 

1 

1 

1 . . . 


Consider the System Below 


'1 

0 

0 

0 

0 

o . . 3 



" o ' 

0 

2x 

0 

0 

0 

0 ... 

A 


X 

x 2 

0 

3x 2 

0 

0 

0 ... 

A 


x 2 

.x 3 

x 1 

0 

4x 3 

0 

0 ... 

A 


X 3 

X 1 

X* 

x J 

0 

5x J 

0 ... 

A 

= 

x J 

x 5 

x 5 

x 5 

x 5 

0 

6x 5 . . . 

A 


X 5 

^ • 





V 

v ■ ^ 




Let F(x) = l+f 1 x+f 2 x"’+f 3 x 3 +... and C' (x) = x+x 2 +x 3 +... be ordinary generating function. 
fi=0 


( 1 ) 


* 


2f 2 x=x 
f!X 2 +3f 3 X 2 =X 2 
f!X 3 +f 2 X 3 +4f 4 X 3 =X 3 
fix 4 +f 2 x 4 +f 3 x 4 +5f 5 x 4 =x 4 , * 

Summing on both sides of (*) we get a first order differential equation 
F' (x) +fi(x 2 +x 3 +x 4 +...) + f 2 (x 3 +x 4 +x 5 +...) +... = x+x 2 +x 3 +x 4 +x 5 +....=C' (x) 

F' (x) +fix(x+x 2 +x 3 +x 4 +...) + f 2 x 2 (x+x 2 +x 3 +x 4 +x 5 +...) + f 3 x 3 (x+x 2 +x J +x 4 +x 5 +...)... = C' (x) 

F' (x) + ff[X + f 2 x 2 + f 3 x 3 + )(x+x 2 +x 3 +x 4 +x 5 +...) = C' (x) 

F (x) + (f , x + f 2 x 2 + f 3 x 3 + ) C' (x) = C' (x) 

F’(x)+(F(x)-l)C’(x)=C’(x) => F’(x)= -(F(x)-l)C’(x)+C’(x)=(l-F(x)+l)C’(x) =(2-F(x)C’(x) = C'(x) Solving 

the differential equation with F(0)= 1 and C(0)=0 

we have F(x ) = 2 — e~ c ** 
Now C’(x) = => C(A) = —x — ln( 1 — x). substituting this in ** we have 

F{x ) — 2 — e x + xe x 
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But we know that e 


X 1 x 3 x" 

1 + X H 1 h ... H h . 

2! 3! n\ 


. X 2 X 3 x" \ X 2 x 3 x n 

F(x) = 2 — \ 1 + X H 1 h...H h l+xfl + xH 1 h...H h ) 

1 2! 3! n\ 2! 3! n\ 





=1 + - 1 )^ n =1 + In-i(kv)*" hence /„ = — 

^ n -lA( n _i)i n ,/ ^n- jn n] 

Now we want to evaluate the upper left comer determinants as follows 


|0| =0, 


10 



0 

l 

0 


0 

1 

0 

0 

ii 

= -l, 

i 

0 

2 

= 2, 

1 

0 

2 

0 

k 

ol 


0 

1 

1 

0 

3 



l 

1 


1 

1 

1 

0 


= -3,. 


where the n th such determinant will be denoted by D n n > 1 
One way to determine say D 5 is as follows, consider the system 


1 0 0 0 0 


0 2 
1 0 
1 1 0 
1 1 1 


0 0 0 
3 0 0 
4 0 
0 


r f \ 


/i 

fi 
/ 3 
f* 

J\fsJ 


1 

1 

1 


where the right hand side is the first column from the original matrix of (1) for n=5 
By crammer’s rule and properties of determinants, we have 


fs ~ 


1 

0 

0 

0 

0 

0 

2 

0 

0 

1 

l 

0 

3 

0 

1 

l 

1 

0 

4 

1 

l 

1 

1 

0 

1 

r 

0 

0 

0 

~0 

0 

2 

0 

0 

0 

1 

0 

3 

0 

0 

1 

1 

0 

4 

0 

1 

1 

1 

0 

5 


: 4 =((— l) n 1 — y.In general by induction we have n!f„ =(-l) n_1 D n 


( 2 ) 


i.e we have D n — (—1)™ 1 n! f n =(—!)" x (n — 1) for all n>0 
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Proposition: suppose 


D 



1 

0 

0 

0 

0 

0 . . . i 

b 2 

b l 

2 

0 

0 

0 

0 . . . 

b 2 

b 2 

b i 

3 

0 

0 

0 . . . 

b 4 

b 2 

b 2 

b > 

4 

0 

0 . . . 

b 5 

b A 

b 2 

b 2 

b i 

5 

0 . . . 

b e 

b 5 

b 4 

b 2 

b 2 

b i 

6 : 

b l 

b e 

b 5 

b 4 

b 2 

b 2 

b x . . . 

b S 

b i 

b 6 

h 5 

b 4 

b 2 

b 2 . . . 

V • 






. y 


n— 1 


n = 1 


Let B' (x) = ^ b n x" 1 and let F(x) =1+ ^ f n X " be the generating function for some number such that F' (x) 
(2-F(x)) B'(x).Then 

_ r (-1 r- L fr 

n If— l) n_1 n! : 


D n 


(— 1 )" 1 f n if F(x)is an EGF 
[(— l) n-1 n!/ n if F(x)is an OGF 


Proof: If 


D = 


(A 

1 

0 

0 

0 

0 

0 . . . i 

b 2 

bl 

2 

0 

0 

0 

0 . . . 

b 2 

b 2 

b, 

3 

0 

0 

0 . . . 

b 4 

b 2 

b i 

b i 

4 

0 

0 . . . 

b 5 

b 4 

b 2 

b 2 

b , 

5 

0 . . . 

be 

b e 

b 4 

be 

b 2 

b i 

6 . . . 

b 7 

b 6 

b e 

b A 

be 

b 2 

b t • • • 

b i 

b 7 

be 

b e 

b 4 

be 

b 2 . . . 

V ■ 






• y 


D j = b !, D 2 = 


1 

b i 


- b| - b 2 , D 3 — 


1 

b i 


b } b 2 b { 


= b 1 3 -3b 1 b 2 +2b 3 ... 


If F(x) is an E.G.F of f x , f 2 then consider the system 


f i 

0 

0 

0 

0 . 

' /, i 


(bA 

b l x 

2x 

0 

0 

0 . . 

fj 2! 


b 2 

b 2 x 2 

b,x 2 

3a 2 

0 

0 . . 

FI 31 


b 2 

b 2 x 3 

b 2 x 3 

b ; x 3 

4x 3 

0 . . 

fju 

= 

b 4 

b 4 x 4 

b 2 x 4 

b 2 x 4 

b x x 4 


fj 5! 


b 5 


/V 


J V • J 
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Then /i + 2 ^-x + 3^-x 2 + — I- AC^i* + b 2 x 2 + •••) + ^-(hiX 2 + ^ 2 * 3 + •••..) = C (x) 
F(x ) + f ± xC (x~) + ^-x 2 C(x) + ^x 3 C'{x) + "-C (x) 

F '(x) + (A* + ~x 2 + ^x 3 )C'(x) + ••■=C'(x) 

F (x) + (F(x) - l)C’(x) = C (x) 

F (x) = (2 — F(x))C (x) 


And by using Cramer’s rule on ***** we have 


fi — b 1 — D lt 


k 

2 ! 


h_ _ 

3! ~~ 


1 0 b ± 

2 &2 
&2 bi b 3 
1 0 o' 

&! 2 0 
b 2 t>i 3 


1 1 

fell 

l&l 

1 b2~bi Z 

1 1 

°i _ 2 

l&l 

2| 


2 

b-y I 

I , |&1 2 


*1 

+£?1 , , 

1 1 |&2 


-6 


=* /2 — 


2 — 3 b ib 2 +&A 

-6 


£3 

3! 


=> f3 — D 3 

Continuing in this way we have D n — (— 1 ) n_1 / n 


Example 1 


Consider 1,2,4,8,16,32,64, = fo, fi, f2> f 3 , ■■■ which has an O.G.F 

F(x) = 1 + 2x + 4x 2 + 8x 3 + 16x 4 + 32x 5 + 64x 6 + 

K J 1-2X 


F'(x) = - — And 2 — F(x) = hence B (x) = 

v 2 (1— 2x) 2 V 7 1-2* V 7 2-F(x) 


(— ) (— ) 
Vl— 2x/ Vl— 4x/ 


=> B\x) = 2 + 12x + 56x 2 + 176x 3 + 992x 4 + ••• 


Hence \ = 2 ,b 2 = 12 ,b 3 = 56 ,h 4 = 176, b 5 = 992,... 


Impact Factor(JCC): 2.7341 - This article can be downloaded from www.impactjournaIs.us 



Masreshaw Walle Abate 


54 


2 

1 

0 

0 

0 

0 

0 . . . 

12 

2 

2 

0 

0 

0 

0 . . . 

56 

12 

2 

3 

0 

0 

0 . . . 

176 

56 

12 

2 

4 

0 

0 . . . 

992 

176 

56 

12 

2 

5 

0 . . . 


992 

176 

56 

12 

2 

6 . . . 



992 

176 

56 

12 

2 .. . 




992 

176 

56 

12 . . . 

v ■ 






. y 


then D n — (—1)™ 1 n\ f n = (—1)" 1 2 n n! for all n > 1 

Example 2: Let F(x) = ^ ^ be an ordinary Generating Function for a certain sequence. 


=> F(a) 


4 

(l-4x) 2 


c'CO = 


F (x) (l-4x) 2 (l-4x) 2 

2-F{x ) ~~ 2— ( — - — ) _ 2 ~ 8 *~ 1 ) 
K 1-4X J 1-4X J 


= (— ) (— ) 

V1-4X/ 


=4(1 + 4a + 16x 2 + 64x 3 4- 256x 4 + 1024x 5 + 4096x 6 + - )(1 + 8x + 64x 2 + ••• ) 
=(4 + 16a: 4- 64a: 2 4- 256a 3 4- 1024a 4 4- 4096a 5 4- - )(1 4- 8a 4- 64a 2 4- 512a 3 4- ■•• ) 
==(4 4- 48a 4- 448a 2 4- 3840a 3 4- ••• by product rule. 

Where b x = 4 ,b 2 = 48, b 3 = 448 ,b 4 = 3840, ... 


Now if 


D = 


' 4 

1 

0 

0 

0 

0 

0 . . . ' 

48 

4 

2 

0 

0 

0 

0 . . . 

448 

48 

4 

3 

0 

0 

0 . . . 

3840 

448 

48 

4 

4 

0 

0 . . . 


3840 

448 

48 

4 

5 

0 . . . 



3840 

448 

48 

4 

6 . . . 




3840 

448 

48 

4 . . . 





3840 

448 

48 . . . 

V 






. ^ 


then D n — (— l) n_1 n! f n = (— l) n_1 4 ”n! for all n > 0 
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